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Abstract 

Measuring the Higgs self-coupling is one of the crucial physics goals at the LHC Run-2 and 
other future colliders. In this work, we attempt to figure out the size of SUSY effects on the 
trilinear self-coupling of the 125 GeV Higgs boson in the MSSM and NMSSM after the LHC Run- 
1. Taking account of current experimental constraints, such as the Higgs data, flavor constraints, 
electroweak precision observables and dark matter detections, we obtain the observations: (1) In 
the MSSM, the ratio of has been tightly constrained by the LHC data, which can be 

only slightly smaller than 1 and minimally reach 97%; (2) In the NMSSM with A < 0.7, a sizable 
reduction of Af^ can occur and minimally reach 10% when the lightest CP-even Higgs 

boson mass mh^ is close to the SM-like Higgs boson m/ij due to the large mixing angle between 
the singlet and doublet Higgs bosons; (3) In the NMSSM with A > 0.7, a large enhancement or 
reduction —1.1 < < 2 can occur, which is accompanied by a sizable change of 

/iir+T“ coupling. The future colliders, such as the HL-LHC and ILC, will have the capacity to 
test these large deviations in the NMSSM. 
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I. INTRODUCTION 


, 3 , 


much effort has been devoted 


With the discovery of the Higgs boson at the LHC [1 
to study its properties. So far, the measurements of its couplings and quantum numbers are 
compatible with the standard model (SM) predictions at l-2a level. However, to ultimately 
understand its nature, we need to fully reconstruct the Higgs potential at the LHC and 
future e’''e“ colliders |3|, 1^. The parameters in the Higgs potential determine the relations 
among the Higgs masses and self-couplings. Measuring these relations is therefore crucial 
for our understanding of the Higgs nature. 

In the SM the tree-level Higgs potential is given by 


^ (0, t> + hf 

which yields the following trilinear and quartic self-couplings 


( 1 ) 
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Here v = ~ 246 GeV is the vacuum expectation value of the Higgs field and 

rrih — 125 GeV is the Higgs boson mass. Within the SM, the trilinear Hi^s coupling receives 
the dominant correction from the top quark loop, ^Af^ ~ mf / 5|, which reduces its 
tree-level value by about 10%. Hence, the determination of the Higgs trilinear coupling Xhhh 
and quartic coupling Xhhhh can directly test the relation in Eq. ([2]) which is obtained from the 
minimization of the Higgs potential. At the LHG, the only way to measure the Higgs trilinear 


coupling is through the Higgs pair proc 
mechanism and has a small cross section 


uction, which is dominated by the gluon fusion 
Gj. However, in many new physics models, such as 
the minimal supersymmetric model (MSSM) and the next-to-minimal supersymmetric model 
(NMSSM), the Higgs pair production rate can be significantly altered by new particles and 
Higgs couplings [7|. Among various decay channels of the Higgs pair, the 46 final state 
has the largest ftaetion g but the tare process M, .ly, g expected to have the most 
promising sensitivity due to the low backgrounds at the LHG [9| . The recent applications of 
jet substructure and other techniques to Higgs pair production have been found to improve 
the sensitivity to the trilinear Higgs couplings in and W^W~ final states [^. On 

the other hand, the measurement of the Higgs quartic coupling is more challenging due to 
a much smaller cross section of triple Higgs production at the LHG. 
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The supersymmetric (SUSY) corrections to the Higgs self-coupling have two kinds of 
sources: one is the mixing between the Higgs bosons, and the other is radiative quantum 
effects. For the first kind, the authors in ll| studied the Higgs self-coupling in some simpli- 


hed SUSY models, while in 


12| the Higgs self-couplin gs i n the MSSM and NMSSM (with 


a decoupled singlet boson) were investigated. Also, in [l3| the authors studied the proper¬ 
ties of the Higgs bosons in the NMSSM with A > 0.7 (called A-SUSY) and found a sizable 
enhancement in the Higgs self-coupling. For the second kind of SUSY corrections, the loop 


corrections to the Higgs self-couplings have 


22| or Feynman diagrammatic approach 


23 


Deen studied using the effective potential flJ - 


25| in the MSSM and NMSSM. Recently, the 


leading two-loop SUSY-QCD corrections from the top/stop sector in the MSSM have been 


performed 


26l |. Since all these previous studies are limited to some simplihed or special 


cases and the relevant experimental constraints are not fully considered, in this work we 
give a comprehensive study for the SM-like Higgs self-coupling in the MSSM and NMSSM 
with both A < 0.7 and A > 0.7 by considering all the relevant experimental constraints after 
the LHC Run-1. 

The existing experimental data, both from low energy precision measurements and high 
energy direct searches, may have imposed important constraints on the Higgs self-coupling 
in SUSY models. For example, in the MSSM, the SM-like Higgs self-coupling is sensitive 


to the pseudo-scalar mass rriA and tan (3, which con 


have been tightly constrained by the 


HC direct search of a light non-SM Higgs boson |27l-l29[| as well as from precision R-physics 
3ol |. Furthermore, the measured mass of the SM-like Higgs boson requires rather heavy 
stops and/or large Higgs-stop trilinear couplings, and the direct searches for the stop pai r 


production have also pushed stop masses above hundreds of GeV in the natural SUSY 




Since the Higgs self-couplings are sensitive to stop masses and Higgs-stop trilinear couplings, 
all these constraints should be taken into account. 

The structure of this paper is organized as follows. In Section m we will briefly describe 
the Higgs sectors of MSSM and NMSSM. In Section UTTl we perform a scan over the parameter 
space of each model and present the numerical results for the trilinear self-coupling of the 
SM-like Higgs. Finally, we draw our conclusions in Section HVl 
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II. HIGGS TRILINEAR SELF-COUPLINGS IN MSSM AND NMSSM 

A. Higgs trilinear self-couplings in the MSSM 

In the MSSM there are two doublets of complex scalar helds with opposite hypercharges: 



The scalar Higgs potential consists of the H-terms and F-terms of the superpotential as well 
as the soft SUSY-breaking mass terms. Among them, the F-terms determine the quartic 
Higgs interactions. The full tree-level Higgs potential is given by 

where is the antisymmetric tensor and 2 = d~^ with iriH^ ^ and /i denoting the 

soft SUSY-breaking masses and the higgsino mass, respectively. The parameters mi ^2 can 
be eliminated by the minimization condition of the Higgs potential, while the parameter 
is traded for the pseudoscalar mass Ma- The quartic Higgs couplings are hxed in terms of 
the SU{2) X U{1) gauge couplings g and g' in the MSSM. 

After the electroweak symmetry breaking, the neutral components of the two Higgs helds 
^ develop vacuum expectation values (vevs) Vu,d and can be decomposed into scalar and 
pseudoscalar components as 

Re = {vd + Hc^ - hso^)/V2, Im - As/3)/V2, (5) 

Rei/o = K Hs^ + hc^)/V2, ImH^ = (G% + Ac^)/V2 (6) 

where h, H and A are the neutral physical Higgs bosons and is the would-be Goldstone 

boson. The vevs are dehned as Vu = vs/s and Vd = vcp with v ~ 246 GeV (here and in the 

following we use the notation = cosx, = sinx). 

Taking the third derivatives of ploYSSM) respect to the physical Higgs helds yields 
the trilinear Higgs couplings. In the physical mass eigenstates, the neutral GP-even Higgs 
trilinear couplings at leading order are given by 

.{0,MSSM) _ 3M| .{0,MSSM) _ Ml 

'^hhh ~ ^ C2q.Sq,_|_^, -^Hhh ~ ^ [^^2aSa+l3 C2qCq_|_^J, 

.(o,MSSM) _ 3M| .(o,MSSM) _ M^ 

'^HHH — ^ C2aCa+l3j '^HHh — ^ [^S2aCa+/3 “r C2aSa+/3\- (.oj 
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In the MSSM, either the lighter scalar h or the heavier scalar H can be the SM-like Higgs 
boson. The latter interpretation occurs for low values of Ma (between 100 and 120 GeV) with 
moderate values of tan/3 (about 10). In this case, H has approximately SM-like properties, 
while the other four Higgs bosons of the MSSM would be rather light and have a mass of 
order 100 GeV or even below. A dedicated scan for this region of parameter space has been 
performed in 32| and it was found that this scenario can be excluded by recasting the LHG 
search for H/A —)■ 27|. In addition, the latest ATLAS limits from searches have 

also excluded such a possibility j^. Thus, in this work, we will only study the case that 
h is the SM-like Higgs boson in the MSSM. According to Appelquist-Garazzone decoupling 
theorem, the MSSM must go back to the SM in the decoupling limit. In the tree-level Higgs 
sector, we can make this limit by setting —)■ cxo, which gives a —)■ /3 — Applying this 

relation to the hrst identity in Eq. ([7]) yields 


A 


iO,MSSM) 


3M| 2 


hhh 


G/3 — 


^7T) ^ 


(9) 

V V 

where — MzC 2 p is the lighter GP-even Higgs mass at tree-level. This demon¬ 

strates that the lighter Higgs boson h in the MSSM almost behaves like the SM Higgs boson 
in the decoupling limit (even when the loop corrections are included) 23l. l24l|. 


B. Higgs trilinear self-couplings in the NMSSM 


After the discovery of a 125 GeV Higgs boson, the NMSSM 33(] seems to be more favored 
than the MSSM because it can naturally give such a Higgs boson without very heavy top- 
squarks 3^. More importantly, this model can solve the /i-problem: after the singlet held 
develops a vev {S) = Vs/\/2, an effective /i-term (/ie// = Ans/\/2) is dynamically generated. 
Due to the contribution of the singlet scalar held S, the full tree-level Higgs potential can 
be written as 


yloMMSSM, ^ (|AS|2 + „j5,J/fX + (|ASr+r4ji/ti/, + m||Sp 
+l(9l + - H\H,f + 

+ + (tSy + + h.c] , (10) 

o 

where k and A are dimensionless parameters, and A\ and are the corresponding trilinear 
soft breaking parameters. To clearly show the properties of the Higgs sector, we can expand 
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33| 


the neutral scalar fields around the vevs as 

Rehr° = (urf — ifsiu/d + hcos/9)/\/2, Imif^ = (Psin/3 + G°cos/9)/\/2, 

ReP° = + Pcos/5 + hsin/9)/\/2, IniP° = (Pcos/9 — G°sin/3)/\/2, 

Re^ = (u,+ s)/\/2, lmS = Ps/^. (11) 


Substituting Eq. flTTl) into Eq. flTOl) . we obtain the mass matrix squared M| for the neutral 
CP-even Higgs bosons and the trilinear Higgs self-interactions as 

/P 




h 

VI 


P ^hahahjhoihfjh^ , 


with ha,j 3 ,'Y = H, h, s. The tree-level Mg- - are given by |35[ | 


( 12 ) 



= Ml + {Ml - ^A^u^) sin^ 2/3, 

(13) 

Ms 12 

= -i(M|-iAV)sin4/3, 

(14) 

M's 13 

= — \/2Au/ix cot 2/3, 

(15) 

Ah's 22 

= M| cos^ 2/3 +-A^u^ sin^ 2/3, 

(16) 

Ah's 23 

= V2Xvfi{l — x), 

(17) 

M's 33 

2 ^4 A^u^ K,X 0 . 

= + ^A^Ia+ 2 X 2 Psm2/3, 

(18) 


with 


Ml = (VAx + KVs) , x=^{Ax + 2V). (19) 

sm2p V / 2/i A 

Here, it should be mentioned that the mass parameter Ma in the NMSSM becomes the mass 
of the heavy pseudoscalar Higgs boson only in the MSSM limit (A, k —?• 0 with the ratio k/X 
fixed). In the NMSSM, Ma can be traded by the soft parameter Ha- 

The CP-even Higgs mass eigenstates hi {i = 1,2,3) can be obtained by diagonalizing M| 
with a rotation matrix O 


K^OicK, {K = H,h,s), diag(mj ,,7ii| ) = 0M|0 


( 20 ) 


with Oia being the elements of the rotation matrix satisfying the sum rules 


I (02 I (92 _ 1 


3a 


( 21 ) 
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The mass eigenstates hi are aligned by the masses mh^ < mh 2 < '^h^- The singlet or non- 
SM donblet components in a physical Higgs boson hi is determined by the rotation matrix 
elements Oi(^H,s)- With Eq. (I2U]) . the corresponding tree-level trilinear Higgs conplings in the 
mass eigenstates hi are given by 


, {0,NMSSM) _ . 

'^hihjhk ~ ^ia^j^^k'y'^hahph^- 


( 22 ) 


In the NMSSM, we take hi or /i 2 as the SM-like Higgs boson when |C>(i, 2 )/iP > 0.5. In 
general, dne to the introdnction of the singlet s and its couplings to the MSSM Higgs sector, 
the mass of the SM-like Higgs boson can be lifted by the extra large A-term at tree 


level, as shown in Eq. (fT 6 |) . The value of A at the weak scale is upper 


order for the NMSSM to remain perturbative up to the GUT scale 


M, 


Dounded by 0.7 in 


35j. Whereas, the 


case of A > 0.7 (dubbed as A-SUSY model) is still of interest because it can suppress the 
sensitivity of the Higgs mass with respect to changes of the soft SUSY-breaking masses and 
keep the hne tuning at a moderate level even for stop masses up to 1 TeV. So, in our study, 
we consider both A < 0.7 and A > 0.7 cases: 

• For A < 0.7, in addition to the tree-level A contribution, the mixture of the singlet 
s with the MSSM Higgs h and 77, as shown in Eq. flT^ . in particular with h, could 
further modify the SM-like Higgs mass. If H is decoupled, when the 

mass eigenvalues for the SM-like Higgs boson is pushed up by the positive mixing 
effect after the diagonalization of the h — s mass matrix. However, when , 

the mass eigenvalues for the SM-like Higgs boson mh^ is pulled down by the negative 
mixing effect. Without the negative mixing effect, the maximal tree-level SM-like 


Higgs mass rrih^ can only reach to about 110 GeV |3g, l37|, which needs a sizable loop 


correction from the stop sector to obtain a 125 GeV Higgs boson like in the MSSM. In 
this sense, h .2 being the observed SM-like Higgs boson may be more natural than hi in 
the NMSSM. So, we choose the next-to-lightest GP-even Higgs boson h 2 as the SM-like 
Higgs boson in the following discussions for the NMSSM with A < 0.7. Note that the 
lightest Higgs boson hi in this case is predominantly singlet-like and its mass can be as 
light as about 20 GeV in our scanned samples. Gonsequently, the SM-like Higgs boson 
/i 2 can decay into a pair of light scalars hi and hence the 77 and ZZ* signal rates are 


suppressed. In order to be consistent with the LHG 
of 7.2 —)■ hihi was found to be less than about 30% 


diggs data, the branching ratio 


38| and may be tested through 
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hihi —)■ bbfi~^n~ production channel at the 14 TeV LHC [39|. However, due to our 
interests in the large mixing region, we only display the results with rrih^ > m/i 2/2 in 
the following analysis. We will also decouple the heaviest CP-even Higgs boson h^, by 
requiring rrih^, > 1 TeV and focus on the singlet-doublet system. On the other hand, as 
mentioned above, if hi is the SM-like Higgs boson, the value of A tends to be large in 
order to maximize the tree-level Higgs mass (so in our study for A-SUSY with A > 0.7, 
we will choose hi as the SM-like Higgs boson). This feature may lead to a sizable 
change in Higgs self-coupling when doublet-singlet mixing effect is large. We checked 
this possibility and found that the ratio can vary from 0.29 to 1.17 for 

our samples with A < 0.7. 

For A > 0.7, we choose the lightest CP even Higgs boson hi as the SM-like Higgs boson. 
The reason is that for A > 0.7 the tree-level Higgs mass will be signihcantly lifted. If 
h 2 is assumed as the SM-like Higgs boson, the large A-term in and the positive 
doublet-singlet mixing effect can readily make the SM-like Higgs mass rrih^ exceed 125 
GeV. Thus, such a choice is strongly disfavored by the LHC observed Higgs mass 0 
and will not be further studied in our work. If hi is the SM-like Higgs boson, as pointed 
before, the cancelation between the tree-level A-term and the negative doublet-singlet 
mixing effect can easily yield a 125 GeV SM-like Higgs boson. However, A can not be 
too large, because the large trilinear and quartic couplings of singlet scalar held s will 
largely contribute to the scattering amplitude ss ss. In the high energy limit, the 
unitarity condition requires |A| < 3 and |k| <3. Furthermore, if combined with the 
dark matter relic abundance constraint, the above-mentioned unitarity bound can set 
a generic upper bound 20 TeV for the heavy Higgs masses 40(]. The requirement of 
perturbativity up to the cut-off scale A, i.e., A(A) < 27r and k(A) < 27r, will set upper 
bounds on A and k at weak scale. In this study, we assume the new unknown strong 
dynamics, for restoring the unitary of scatte ring processes, appears at some scale A 
above 10 TeV, and require A^ -|- < 4.2 


35 




Before closing this section, we note that since both Higgs boson masses and Higgs self¬ 
interactions arise from the Higgs potential, one has to adopt the same method to calculate 
them up to the same order for comparing the self-couplings in SUSY and SM. Although the 
most accurate evaluation of Higgs mass is up to three-loop level in the MSSM there is 
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no corresponding result for the Higgs self-couplings. Hence, in our study, we derive the Higgs 




mass and self-couplings by following the effective potential approaches in the MSSM 
and NMSSM 2l|, where the explicit dominant one- and two-loop corrections to the effective 


potential are presented. Then we coded those expressions in our numerical calculation. On 
the other hand, since the tree-level mixing effects in the Higgs sector are usually dominant 
over the loop corrections in the NMSSM, we will focus on the doublet-singlet mixing effects 
for the NMSSM results. 

III. NUMERICAL CALCULATIONS AND RESULTS 


A. Scan over the parameter space 


In our numerical calculations, we take the input parameters of the SM as 431 


mt = 173.5 GeV, mw = 80.385 GeV, mz = 91.1876 GeV, 

= 4.18 GeV, a^i^Mz) = 0.1184, a{mz)-^ = 128.962 (23) 


We use NMSSMTools-4.4.1 


44j | to perform a random scan over the parameter space. Note 


that for any given value of p-e//) fhe phenomenology of the NMSSM is identical to the MSSM 
in the limit A, k —?• 0 with the ratio n/ A hxed should be negative and satisfy | A^l < d/tp/A 
to guarantee the squared mass of the singlet scalar to be positive) j^. In our scan for the 
MSSM, we take X = k = 10“^ and = —10 GeV. The validity of such a method has been 


justihed by the authors of the NMSSMTools 


44| and our previous calculations [3^. We also 


numerically checked our MSSM results of nih by using the codes FeynHiggs |^, SOFTSUSY 


46l | and SuSpect and found the results to agree with that given by the NMSSMTools 


within about 1% level when m/j ~ 125 GeV. So it is feasible to use the package NMSSMTools 
in the MSSM limit to study the phenomenology of the MSSM. For simplicity, we decouple 
the sleptons, gauginos and the hrst two generations of squarks by hxing the corresponding 
soft mass parameters at 2 TeV. We also set Mq^ = = Mf /3 and At = Ai, for the third 


generation of squarks. The lower 
by the recent LHG Higgs results 
following: 


imit of tan [3 in the MSSM is taken as 5, which is inspired 


32j. The parameter ranges in our scan are chosen as the 
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(a) For the MSSM, 


5 < taii/S < 60, 0.2 TeV < Ma < 1 TeV, |/x| < 1 TeV, 

0.1 TeV < Mq 3 < 2.5 TeV, < SMq^. (24) 

(b) For the NMSSM with A < 0.7, 

0 < A < 0.7, |fi:| < 0.7, 0.2 TeV < Aa < 1 TeV, 

I^kI < 1 TeV, 1 < tan/3 < 20, 0.1 TeV < /i < 1 TeV, 

0.1 TeV < Mq 3 < 1 TeV, |A| < SMq^. (25) 

(c) For the NMSSM with A > 0.7, 

0.7 < A < 2, 0 < K < 2, 0.2 TeV < Aa < 1 TeV, 

|^k| < 1 TeV, 1 < tan/3 < 20, 0.1 TeV < /i < 1 TeV, 

0.1 TeV < Mq 3 < 1 TeV, \At\ < SMq^. (26) 


In our scan, we consider the following constraints: 


(1) We require the SM-like Higgs mass in the range of 123-127 GeV and consider the 
exclusion limits (at the 95% conhdence level) from LEP, Tevatron and LHC in Higgs 


searches with HiggsBounds-4.2.0 48(]. We also perform the Higgs data ht by calculating 
of the Higgs couplings with the public package HiggsSignals-1.3.0 [49| and require 
our samples to be consistent with Higgs data at 2a level. We choose the SLHA 
input choice of HiggsBounds/HiggsSignals, where the effective Higgs couplings are 
only used to calculate the Higgs production cross section ratios. The Higgs decay 
branching ratios are taken directly from the corresponding decay blocks in the SLHA 
hie generated by the NMSSMTools. 


(2) We require one-loop SUSY predictions for H-physics observables to satisfy the 2a 
bounds as encoded in NMSSMTools, which include B —>■ Bg —)■ B^ 

and B^ Theoretical uncertainties in B-physics observables are taken 

into account as implemented in NMSSMTools. 


(3) We require the one-loop SUSY predictions for the precision electroweak observables 
such as Pi, sin^^gjj, mw and |50| to be within the 2a ranges of the experimental 
values. 
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(4) We require the thermal relic density of the lightest neutralino (as the dark matter can¬ 
didate) is below the 2a upper bound of the Planck value and the spin-independent 
neutralino-proton scatte ring cross section satisfy the direct detection bound from LUX 


at 90% conhdence level 


m- 


(5) We also consider the theoretical constraints from the stability of the Higgs potential 
as encoded in the NMSSMTools. 

B. Results for the MSSM 
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FIG. 1: Excluded region in the tan/? versus plane of the MSSM. 


Precision measurements of the Higgs boson properties (its mass and couplings to other 
particles) at the LHC provide relevant constraints on possible weak-scale extensions of the 
SM. In the usual context of the MSSM, these constraints suggest that all the additional 
non-SM-like Higgs bosons should be heavy. In Fig. [1], we show that in the decoupling limit 
of the MSSM, the region with rriA < 330 GeV has been excluded by various constraints. 
Similar result has been recently pointed out in The lower part of Fig. [1] shows that 
the value of tan/? can not be too small due to the requirement of 125-127 GeV Higgs mass 
in our scan. For a small tan/?, heavy stops (> 10 TeV) or a large mixing parameter At 
is needed to produce a large positive correction to the Higgs mass, which, however, will 
easily lead to vacuum instability and large uncertainty in the numerical calculation. So, we 
focus on irtt^ < 2.5 TeV region in our calculations. The additional constraints on the MSSM 
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parameters imposed by the Higgs data are obtained by using the HiggsBounds-4.2.0 package. 
The key algorithm of HiggsBounds can be described in two steps. Firstly, the HiggsBounds 
uses the expected experimental limits from LEP, Tevatron and LHC to determine which 
decay channel has the highest statistical sensitivity. Secondly, for this particular channel, 
the theory prediction is compared to the observed experimental limits to conclude whether 
this sample is allowed or excluded at 95% CL. With the HiggsBounds, we hnd that most 
of samples (black bullets) with 180 < Ma ^ 330 GeV and 5 < tan /? < 10 have been 
excluded. Particularly, the latest CMS result of searching for H/A —>■ has excluded 

most parameter space with a low Ma and low to moderate values of tan/?, as shown in the 
upper left corner of Fig. [TJ We have checked that the current low energy constraints from 
Bg —)■ and Bg —)■ Xg'y are weaker than H/A —)■ for our interested region (low to 

moderate values of tan/3). We note that the supersymmetric loop corrections generally lead 
to a contribution of the order of a few percent of the SM value. Hence, the suppressions in 
Higgs signal strengthes fi.y.y and /iw* are mostly governed by the increase of the width of 
the lightest CP-even Higgs decay into bottom quarks and tan leptons at low values of tua- 
Fig. [T] suggests that all the additional non-SM-like Higgs bosons should be heavy, with 
masses larger than about 330 GeV. This is the commonly discussed decoupling limit of 
the MSSM. However, as discussed in [5^, it is also possible to have the MSSM parameter 
conditions for alignment independent of decoupling, where the lightest GP-even Higgs boson 
has SM-like tree-level couplings to fermions and gauge bosons, independently of the non¬ 
standard Higgs boson masses. In the alignment region, sin(a — /3) ~ 1 and the bounds on 
the heavy Higgs bosons that arise from the measurements of h —)■ VV may be relaxed. Such 
alignment conditions are associated with very SM-like hti coup ling and tend to be restricted 


to values of tan /3 of order 10 or larger within the MSSM 


54 


56| . As will be shown below. 


for such a large value of tan/3, the supersymmetric contributions to the SM-like Higgs self¬ 
coupling are found to be small (at a percent level). Hence, even in the alignment condition, 
we cannot expect a large deviation of the Higgs self-coupling with the SM value. 

In Fig. |2]we project the samples surviving all the experimental constraints on the planes 
of tua and mi^ versus tan/3. The ratio is always smaller than 1 because of 

the negative MSSM corrections to Higgs self-coupling 23(]. From the left panel it can be 
seen that in most part of the allowed parameter space (the blue triangles), the values of 
AMSsat/asm hhh) are larger than 0.99. When both tan/3 and 
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(TeV) /Wj; (TeV) 


FIG. 2: Scatter plots of the samples surviving all the experimental constraints, projected on the 
planes of tan j3 versus niA and . 


become small, gets smaller, which can minimally reach abont 0.97 for onr 

samples. The reasons for being so close to Af^f^ are the following: (1) The dominant 

MSSM contribntions to X^ff^^/X§^ are from the stop loops, while the Higgs conpling with 
the stops is proportional to 1/sin/?. So it leads to an overall enhancement factor 1/sin^/? 
in the corrections when tan /? is small. However, snch a region is obvionsly not favored by 
the measnred Higgs mass, which needs a large tan /? to enhance the Higgs mass; (2) A light 
rriA canses a large mixing between two CP-even Higgs bosons and can sizably change the 
Higgs conplings with the SM fermions. Bnt as mentioned before, rriA shonld be heavier than 
abont 330 GeV to satisfy the experimental constraints in onr scan. For the right panel it 
shonld be mentioned that since the small valnes of /Xf^ occnr in the small tan/? 

region, heavy stops are nsnally needed to enhance the Higgs mass throngh loop corrections. 

In Fig. [3l we show the MSSM Higgs conplings in comparison with the SM predictions. 
The ILC (1 TeV, 1 ab“^) sensitivities to the alteration of the conplings are also ploted, 
where the regions between the bars give too small alterations to be detectable at ILC. The 
HL-LHC (14 TeV, 3 ab“^) sensitivities are mnch worse than ILC and are not shown here. 
From Fig. [3l we have the following observations: (1) In the MSSM, the Higgs gange conplings 
and top-Higgs conplings are respectively changed by the factors sin(/? — a) and cos a/sin 
bnt the ratios Cj^^y^/C^yy and for onr snrvived samples are very close to 

nnity due to uia 3> mz- So, even if these couplings can be measured at percent level at ILC, 
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FIG. 3: Same as FiglJl but showing the Higgs couplings. The ILC (1 TeV, 1 ab sensitivities 


to the alteration of the couplings 


57l | are also ploted (the regions between the bars give too small 


alterations to be detectable at ILC). 


it can not constrain the MSSM parameter space with a small (2) The ratio 

always smaller than one for our samples because a large mixing between the 
stops or heavy stops needed by the Higgs mass interferes destructively with the top loop and 
is suppressed, while can be greater than one due to the constructive 

contribution with the W loop. Although both couplings only slightly deviate from the SM 
predictions, the high precision measurements on at the ILC will be able to exclude 

some part of the parameter space with > 0.977. 


C. Results for the NMSSM (A < 0.7) 

In Fig. m we project the NMSSM (A < 0.7) samples allowed by the constraints (l)-(5) on 
the planes of A versus k and k versus the singlet component \02s\ in the SM-like Higgs boson 
^ 2 . Here we take > mh 2 so that h 2 will not decay into a pair of hi. We can see that 
when A and k approach to zero, tan jS has to be large in order to enhance the Higgs mass. 
While in small tan ^ 10 region, the values of A for most samples are larger (in magnitude) 
than K. Since the singlet component \02s\ in the SM-like Higgs boson h 2 can potentially 
affect the Higgs self-coupling in Eq. fl22l) . we also show the dependence of \02s\- Since the 
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FIG. 4: The NMSSM (A < 0.7) samples surviving all the experimental constraints, projected on 
the planes of A versus k and the singlet component \02s\ in the SM-like Higgs boson /i 2 . 


singlet component \02s\ in the SM-like Higgs boson h 2 can potentially affect the Higgs self¬ 
coupling in Eq. fl^ . we also show the dependence of \02s\ on A and k for the NMSSM (with 
A < 0.7) samples which survive all the experimental constraints. It can be seen that \02s\ 
deceases if both A and k go to zero. However, if A 3> |k| ~ 0, which usually happens for 
small tau/d < 10, \02s\ can be sizeable and lead to a large deviation from the SM prediction 
in triple Higgs boson coupling. For example, for A = 0.017 and k = —0.0038, O 2 S ~ 0.215 
and 0.7. The reason is that in the singlet-doublet system, after decoupling 

hs by requiring Mh^ > 1 TeV, \02s\ is proportional to sin0, where the mixing angle 9 that 
determines the mixture of the singlet and the SM-like Higgs states can be approximately 
expressed as tan 20 ~ 2M|3/(M|2 — Mig) 35, 37|. Here, are the CP-even Higgs mass 
matrix elements listed in Eqs. fll3tll8p . In case that A is not too small, a large 9 can happen 


when M 22 — Tfig ~ Mig. This leads to a large singlet component in the SM-like Higgs state 
and hence a sizeable modihcation of the Higgs trilinear couplings. 

In Fig. [5], we display the dependence of /X§^ versus A and tan/? for the NMSSM 

(A < 0.7). As mentioned above, due to the Higgs mass constraint, most of the allowed 
model samples tend to have large values of A. Furthermore, the ratio /X^jf is not 

sensitive to the value of tan/3. We also show the dependence of /Xf^ versus the 

singlet component \02s\ in the SM-like Higgs boson h 2 and the lightest CP-even Higgs mass 
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FIG. 5: Same as Fig. 01 but showing the dependence of versus A and the singlet 

component \02s\ in ^ 2 ; where m/jj < 2m/ij. The ILC(1 TeV, 1 ab“^) and HL-LHC(14 TeV, 3 ab“^) 
sensitivities are also plotted (the region below each horizontal line is detectable). 


nihi for fho NMSSM (A < 0.7). Besides, we plot the expected ILC(1 TeV, 1 ab“^) and HL- 
LHC(14 TeV, 3 ab“^) sensitivities to Higgs self-conpling from the direct measnrements j^. 
We can see that /^Ih becomes small with the increase of the singlet component 

\02s\ and can minimally reach 0.1 in the allowed parameter space. Meanwhile, snch a large 
mixing can make the mass of the lightest singlet-dominant CP-even Higgs boson hi close to 
the 125 GeV SM-like Higgs boson ^ 2 . 

In Table [H we present the properties of two CP-even Higgs bosons hi and h 2 for snch 
a benchmark point at 14 TeV LHC. We can see that the cross sections of the single hi 
prodnction are close to those of SM-like h 2 becanse of the large doublet and singlet mixing 
components in both hi and h 2 - However, the branching ratio hi VViV = Z, IVy,^') 
is greatly suppressed by the increase of the partial width of hi —)■ bb. Thus, the observed 
production rate of gg ^ hi ^ VV is much smaller than that of SM-like h 2 . We also checked 
and found that although the cross section of gg ^ hi ^ can reach 2.06 pb, it is 

still smaller than the upper bound given by the current LHC searches for H/A ^ 

28| . Besides, it should be mentioned that the sizable modihcation of the self-coupling of the 
SM-like Higgs boson h 2 always accompanies with the great changes in other Higgs couplings, 
which can be seen in Fig. 01 So, given the limited sensitivity of measuring the Higgs boson 
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TABLE I: A benchmark point for the NMSSM with A < 0.7 (/i 2 is the SM-like Higgs boson). The 
cross sections (pb) are calculated for LHC-14 TeV. 


A 

K 

tan/3 

p (GeV) 

AIa (GeV) 

(GeV) 

\NMSSM/\SM 
^3/12 / ^Zh2 

0.16 

0.31 

17.9 

105.3 

1461.9 

-716.2 

0.798 

ruhi (GeV) 

^gghi 

O'VVhi 

O'Whi 

crzhi 

^tthi 

^bbhi 

119.56 

25.41 

2.14 

0.85 

0.49 

0.33 

0.54 



Bvh^^zz* 

Brhi^WW* 


B^hi-ybb 


0.135% 

3.60% 

0.717% 

7.88% 

2.12% 

77.3% 

8.10% 

mh 2 (GeV) 

^ggh2 

0'VVh2 

0'Wh2 

(^Zh2 

^tth2 

^bbh2 

127.30 

24.81 

2.05 

0.71 

0.41 

0.29 

0.14 

>77 

Brh2^gg 

Bvh^^zz* 

Brh2^WW* 

^'^h2—>CC 

BT'h2^bb 


0.381% 

8.11% 

3.927o 

34.0% 

4.14% 

44.6% 

4.55% 


self-couplings, we anticipate that the precision measurement of Higgs couplings to gauge 
bosons and fermions could test this scenario at the LHC and future colliders prior to the 
direct detection of triple or quartic Higgs couplings. 



FIG. 6; Same as Fig. [3]but for the NMSSM (A < 0.7). 

In Fig. [6l we show the Higgs couplings in the NMSSM (A < 0.7) surviving all the exper- 
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imental constraints. For comparison, we also show the discovery potential of the expected 
ILC(1 TeV, 1 ab“^) and HL-LHC(14 TeV, 3 ab“^) ^|. From this hgnre, we obtain the 
following observations: (1) Dne to the singlet admixtnre in the SM-like Higgs boson, both 
Higgs gange conplings and top-Higgs conpling can be maximally rednced by abont 30% in 
the allowed region, which is mnch larger than in the MSSM. So, the expected measurements 
of the Higgs gauge couplings at the HL-LHC and ILC can exclude the parameter space 
with < 0.82 and < 0.93, respectively; (2) With the increase 

of the singlet component in the Higgs couplings, both and 

are signihcantly reduced. On the other hand, due to the additional tree-level contribution 
(~ An sin 2/?) and the positive mixing effect, we hnd that a stop with mass less than 200 GeV 
is still allowed by the SM-like Higgs mass constraint in the NMSSM (similar results have 


been obtained in previous NMSSM works ^]). Consequently, the ratio be¬ 

comes larger than one, due to the constructive contribution from the light stop in loop. We 
also note that even when approaches to can still be enhanced 

by about 


D. Results for the NMSSM (A > 0.7) 


In Fig. we display the dependence of versus tan/3 and for the 

NMSSM (A > 0.7). Similar to Fig. HJ the larger A is, the smaller tan/3 becomes to satisfy 
the requirement of the Higgs mass. The ratio can vary from -1.1 to 1.9 in our 

scan. For example, the ratio /Xfff is equal to 1.89 when A = 1.51 and k = 0.67, 

and is —1.04 for A = 1.57 and k = 1.16. Since we require irth^ > 1 TeV, for our samples with 

is approximately proportional to A^ and becomes large with the 
increase of A. In our scan ranges, such a feature could either enhance or suppress the Higgs 


self-coupling with respect to t 
pair production cross sections 


re SM prediction, and yield potentially large effects in Higgs 
13 


■ l59[| . However, in general, the two masses mh 2 and rrih^ are 
virtually independent and the mixing patterns are complicated. It is worth mentioning that 
the value of ^ 1 strongly relies on the mass of the next-to-lightest CP-even 

Higgs boson h 2 . To be specihc, when n becomes small (large), rrih^ is inclined to be light 
(heavy). If rrih^ is lighter (heavier) than about 400 GeV, /Xf^ is smaller (larger) 

than unity for most samples. The properties of h 2 will be discussed in the following. We 
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0.74 

0.94 

1.1 

1.3 

1.5 

1.7 
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FIG. 7: The NMSSM (A > 0.7) samples surviving all the experimental constraints, showing the 
dependence of versus tan/3 and m/ij- The ILC(1 TeV, 1 ab“^) and HL-LHC(14 

TeV, 3 ab“^) sensitivities are also plotted. 


should mention that the large A and k that produce the large deviation of 
jeopardize the perturbativity up to GUT scale. Thus, the new unknown strong dynamics 
will appear at some cut-off scale A. On the other hand, the constructions of high scale 
theory by adding vector-like matter can allow for a larger A value and relax the cutoff scale 


60| . which is however beyond the scope of our study. 


to high values 

In Fig. [HI we plot the couplings of ^2 with gauge bosons and r leptons, normalized to 
the SM values. We can see that the gauge coupling h 2 VV is always suppressed due to 
the presence of singlet (s) and non-SM doublet (if) components in h 2 - If ^2 is singlet¬ 
like, the h 2 T^T~ coupling is suppressed as well, while if ^2 is non-SM doublet-like, h 2 r+r“ 
coupling can be enhanced by tan/3. The detailed mixing patterns of s and H in h 2 and 
its couplings have been thoroughly investigated in |^. We checked that the cross section 
gg ^ h 2 ^ T^T~ for our samples is at least one order lower than the current direct search 
bound on the non-SM Higgs bosons 27|, l28|. The main reason is that for a heavy Higgs 
boson ^ 2 , the new decay channels, such as h 2 —>■ hihi 6l|, can be opened and the branching 
ratio of ^2 —)■ will be highly suppressed. 

In Table [TTl we present the properties of two CP-even Higgs bosons hi and /12 for a 
benchmark point at 8 TeV LHC. We can see that all the SM branching ratios of ^2 are 
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FIG. 8: The NMSSM (A > 0.7) samples surviving all the experimental constraints, showing the 
dependence of the Higgs couplings versus m/ij • 


reduced due to the opened decay mode ^2 —)■ hihi, which can reach 53.6% for mh^ = 282.0 
GeV. Such a Higgs-to-Higgs decay will lead to a resonant SM-like di-Higgs bosons production 
pp —)■ ^2 —)• hihi at the LHC. A resonance feature in the hihi invariant mass can be served 
as a smoking gun to search for the heavy Higgs boson ^ 2 . With one SM-like Higgs boson 
hi decaying to two photons and the other decaying to 6 -quarks, the resonant signal ma y b e 
observable above the di-Higgs continuum background for mh^ < 1 TeV at the HL-LHC 62|. 
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FIG. 9: Same as Fig. [8l but s^ 
TeV, 1 ab“^) sensitivities 57| 


rowing the couplings of 6i. The HL-LHG(14 TeV, 3 ab and ILC(1 
are also plotted. 
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TABLE II: A benchmark point for the NMSSM with A > 0.7 (hi is the SM-like Higgs boson). The 
cross sections (pb) are calculated for LHC-14 TeV. 


A 

K 

tan j3 

p (GeV) 

Ha (GeV) 

H. (GeV) 

\ NMSSM /\SM 
^3h2 /^3h2 

1.54 

0.75 

3.06 

474.00 

1035.41 

-644.30 

-0.203 

nihi (GeV) 

O'gghi 

O'VVhi 


O'Zhi 

^tthi 

^bbhi 

126.9 

42.86 

3.42 

1.17 

0.70 

0.52 

0.29 


Brh^^gg 


Bvh^^ww* 


B^hi-ybb 

Br^^^^+r- 

0.427% 

9.87% 

4.18% 

36.5% 

4.98% 

39.7% 

4.03% 

mh 2 (GeV) 

0'ggh2 

0'VVh2 

(^Wh2 

0 'Zh2 

^tth2 

^bbh2 

282.0 

6.33 

0.90 

0.06 

0.03 

0.03 

0.04 

'®^/i2^77 

h^^h2^gg 

Bvh^^zz* 

Brh2^WW* 

Brh2^bb 

Brh^,^+r- 

^'^h2^hihi 

0.000663% 

0.0173% 

14.1% 

32.1% 

0.0996% 

0.012% 

53.6% 


Next, we present in Fig. |9]the couplings of the SM-like hi to weak gauge bosons and tan 
pair for the NMSSM (A > 0.7). To show their correlation with the Higgs self-coupling, we 
use the red color to highlight the points that satisfy — 1| > 0.3. It can be 

seen that the large Higgs self-couplings corrections correspond to the sizable shifts in Higgs 
couplings with gauge bosons and tan pair. Since the tree-level mass of the SM-like Higgs 
boson hi could easily exceed 125 GeV, hi is likely to have non-negligible singlet and/or 
non-SM doublet components, which makes its couplings deviate from the SM predictions. 
From this hgure, it can be seen that /C^yy is always less than unity since any 

singlet and/or non-SM doublet components in hi will make the Higgs couplings to weak 
gauge bosons smaller than the SM predictions. However, as mentioned above, if the next- 
to-dominant component in hi is the non-SM doublet, the Higgs coupling with the down- 
type fermions may be enhanced and larger than the SM predictions, such as the hiT^T~ 
coupling. Therefore, the future measurements of Ch^r+r- and Ch^vv couplings can give 
strong constraints on the parameters space of the NMSSM with A > 0.7 and set limits on 
the Higgs self-coupling 

From the above discussions, we can see that the large deviation of the Higgs self-coupling 
X^hi for the NMSSM (A > 0.7) is always accompanied by other collider signatures, such 
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as a shift in the Higgs couplings and the production on resonance of the non-SM doublet 
h 2 - At the LHC, the resonant production of ^2 may be observed through the channels 
gg ^ h 2 ^ T~^T~, hihi- However, the sensitivities of these channels strongly depend on the 
mass of h 2 and on its singlet and non-SM doublet components. If ^2 is dominantly singlet, 
both direct searches will not be powerful in probing our scenario at the LHC since all the 
h 2 couplings to SM particles are greatly reduced. Thus, the precision measurement of the 
observed 125 GeV Higgs boson couplings will play an unique role in probing this scenario 
at future colliders. If h 2 is dominantly doublet (non-SM), and if mh^ < mh 2 < 2m/i^, the 
process 5'fi' —t ^2 —>■ r+r“ is still greatly suppressed due to the reduction of the coupling h 2 ti 
despite the fact that the coupling h 2 T^T~ can be maximally enhanced by a factor of 2. On 
the other hand for 2m/i^ < m/jj, the decay ^2 —)■ hihi is open and contributes to the cross 
section of pp —)■ hi hi. Therefore, besides the Higgs coupling measurements, a resonance 
feature in the hi hi invariant mass or an excess in the inclusive hi hi production can be used 
to probe this model at the future LHC. 

IV. CONCLUSION 

We examined the currently allowed values of trilinear self-couplings of the SM-like 125 
GeV Higgs boson (h) in the MSSM and NMSSM after the LHC Run-1. Considering all 
the relevant experimental constraints, such as the Higgs data, the flavor constraints, the 
electroweak precision observables as well as the dark matter detections, we performed a 
scan over the parameter space of each model and obtained the following observations: 

• In the MSSM, the Higgs self-coupling is suppressed relative to the SM value. Such a 

suppression was found to be rather weak and the ratio is above 0.97 due 

to the tightly constrained parameter space, cf. Figs. [Hand [21 

• In the NMSSM with A < 0.7, we consider the case that the SM-like Higgs boson mass 

mh 2 is less than twice of so that ^2 will not decay into a pair of hi. We found that 
the Higgs self-coupling was found to be likely suppressed and the ratio / A™;^ 

can be as low as 0.1 due to the large mixing between singlet and doublet Higgs bosons, 
cf. Fig. [5l In that case, the coupling of the SM-like Higgs boson (^ 2 ) to to W and 
Z bosons and top quark pairs are all suppressed as compared to the SM prediction. 


22 


On the other hand, its couplings to loop-induced processes, such as photon pairs or 
gluon pairs, can be enhanced, cf. Fig. [6l Given the limited sensitivity of measuring 
the Higgs boson self-couplings, we anticipate that the precision measurement of Higgs 
couplings to gauge bosons and fermions could test this scenario at the LHC and future 
colliders prior to the direct detection of triple or quartic Higgs couplings; 

• In the NMSSM with A > 0.7 (also called A-SUSY), the Higgs self-coupling can be 
greatly suppressed or enhanced relative to the SM value (the ratio A™/i can 

vary from -1.1 to 1.9), cf. Fig. 0 when hi is taken as the SM-like Higgs boson. Its 
coupling to W and Z bosons are always suppressed as compared to the SM predictions. 
On the contrary, its couplings to tan pairs can be either enhanced or suppressed relative 
to the SM value, cf. Fig.|9l When mh 2 > it is possible to observe a new resonance 

production in the di-Higgs boson channel at the LHC. While the couplings of ^2 to W 
and Z bosons are always suppressed, its coupling to tan pair can deviate largely from 
the SM value, depending on the mass of ^ 2 , cf. Fig. [HI 

Since the NMSSM can give rather different values (compared with the SM) for the trilin- 
ear self-couplings of the SM-like Higgs boson, the future collider experiments like the high 
luminosity LHC or ILC can probe NMSSM through measuring the Higgs self-couplings. 
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